Reduced fidelity approach for quantum phase transitions in spin-1/2 dimerized 

Heisenberg chains 
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We use reduced fidelity approacli to cliaracterize quantum pliase transitions in tlie one-dimensional 
spin-1/2 dimerized Heisenberg drain in the antiferromagnetic case. Tiie reduced fidelity suscepti- 
bilities between two nearest-neighboring spin pairs are considered. We find that they are directly 
related to the square of the second derivative of the ground-state energy. This enables us to conclude 
that the former might be a more effective indicator of the second-order quantum phase transitions 
than the latter. Two further exemplifications are given to confirm the conclusion is available for 
a broad class of systems with SU(2) and translation symmetries. Moreover, a general connection 
between reduced fidelity susceptibility and quantum phase transitions is illustrated. 

PACS numbers: 75.10.Pq, 03.67.-a, 75.10.Jm 



INTRODUCTION 



Quantum phase transitions (QPTs) is an essential phe- 
nomenon in quantum many-body correlated system. It 
is induced by the ground-state (GS) transition driven 
by external parameters at zero temperature. How to 
characterize QPTs has attracted widespread attention. 
Conventionally, QPTs are described in terms of order 
parameter and symmetry breaking within the Landau- 
Ginzburg paradigm One object in these traditional 
ways is that there is no general method to find the or- 
der parameter for a common system. To overcome this 
problem, a concept called fidelity ^, '3] (see Eqs. ^ and 
([7])) is borrowed from the field of quantum- information 
theory, since it well describes the overlap between two 
states in different phases with different structural prop- 
erties, i.e., order parameters. Thus it dose not need a pri- 
ori knowledge of the order parameter in detecting QPTs. 
It is a purely Hilbert-space geometrical quantity. On 
the other hand, fidelity susceptibility [3, @ (see Eq. ([S])) 
is found more convenient than fidelity itself for its in- 
dependence of the slightly changed external parameters. 
Hitherto, these two connected concepts have succeeded 
in identifying the QPTs of many systems, such as XY 
spin chains and Dicke model XXZ chain ^, Hub- 
bard model 0, Q , frustrated Heisenberg chain 1 1 , Kitaev 
honeycomb model fl^ , extended Harper model ll| . The 
intrinsic relation between the GS fidelity (or fidelity sus- 
ceptibility) and the characterization of a quantum phase 
transition has been unveiled in Ref. p^ . It was shown 
that the singularity and scaling behavior of the GS fi- 
delity (or fidelity susceptibility) are directly related to 
its corresponding derivative of GS energy, which charac- 
terizes the QPTs conventionally. Moreover, the fidelity 
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susceptibility is associated with dynamic structure factor 
for QPTs and with specific heat and magnetic suscepti- 
bility for thermal phase transitions [T3| . 

The above works are all concerned with the global GS 
fidelity. Then there is a natural question that whether 
the fidelity of the subsystem, i.e., the reduced fidelity (or 
named partial-state fidelity) could reflect the QPTs. Re- 
cently, some works have been devoted to this subject. 
Zhou et al. |14i] found that it succeeds in capturing non- 
trivial information along renormalization group flows and 
in detecting the QPTs in XY model [Tsl. Paunkovic et 
al. [16,] showed that it enables them to identify the on- 
site magnetization as the order parameter for the phase 
transition in the conventional BCS superconductor with 
an inserted magnetic impurity system. Kwok et al. [iTj 
tested its effectiveness in charactering the QPTs of the 
isotropic Lipkin- Meshkov-Glick model and the antifer- 
romagnetic one-dimensional Heisenberg model. Mean- 
while, we derived a general expression for the two-site 
reduced fidehty susceptibility (RFS). It has been appUed 
to the study of the Lipkin-Meshkov-Glick model [1^ and 
transverse field Ising model fiol] . We found that the RFS 
shows similar scaling behavior to the global fidelity sus- 
ceptibility. All the above works illustrate that the re- 
duced fidelity approach is also an effective tool in iden- 
tifying QPTs. However, a general relation between RFS 
and QPTs is not established. 

In this work, we apply the reduced fidelity approach 
to the one-dimensional (ID) spin-1/2 Heisenberg chain 
with dimerization. Thanks to the SU(2) and translation 
symmetries, we derive general expressions of the two-site 
RFSs (the RFSs we mentioned below are all for two-site.), 
which are connected closely to the square of the second 
derivative of the GS energy. This result indicates that the 
RFS is an more effective tool to identify the second-order 
QPTs than the second derivative of the GS energy. To 
further testify our conclusion, we exemplify the mixed- 
spin dimerized Heisenberg chain and the spin-1 bilinear- 
biquadratic model as well. These two models are both 
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of SU(2) and translation symmetries too. Furthermore, 
it is illustrated that, in general, the origin for RFS to 
signal QPTs may root in the relation between the reduced 
density matrix (RDM) and the derivatives of GS energy. 

This paper is organized as follows. In Sec. II, we derive 
a general expression of RFS for two Hermitian and semi- 
positive definite density matrices, which are commute 
with each other, and give a direct connection between 
RFSs and QPTs in the dimerized model. In Sec. Ill, 
the critical behavior of the system is studied for both 
finite-size and infinite-size situations. In Sec. IV, two 
further models are enumerated and a possible origin be- 
tween the relation of general reduced fidelity and QPTs 
is illustrated. Finally, a summary is presented in Sec. V. 



II. REDUCED FIDELITY SUSCEPTIBILITY 
AND ITS CONNECTION TO QUANTUM PHASE 
TRANSITIONS 

The dimerized Heisenberg chain is a fundamental spin- 
correlated model. It is of special interest both in the- 
ory and experiment, since it gives a reasonably accurate 
description of many quasi-lD antiferromagnets which 
have two important but structurally inequivalent su- 
perexchange paths that are spatially linked, such as the 
materials of Cu(N03)2 •2.5H2O, (VO)2P207 and various 
aromatic free-radical compounds [2^. Therefore, many 
efforts have been devoted to study its quantum critical 
behavior of the dimerized Heisenberg model using var- 
ious methods, e.g., continuous unitary transformations 
pll |. density matrix renormalization group p^ . concur- 
rence [2^ and block entanglement [l^. Here we employ 
the reduced fidelity approach to study the QPTs of this 
model. 

The Hamiltonian for antiferromagnetic Heisenberg 
chain (AHC) with dimcrization reads 



with 



N/2 



Hr 



— (S2'i-1 • + aS2i ■ S2i+i) , 



(1) 



i=l 



where Si denotes the i-th spin-1/2 operator, and a > 
is the ratio between the two kinds nearest-neighboring 
(NN) couplings. The total number of spins N is required 
to be even and the periodic boundary condition Si = 
Stv+i is assumed. 



A. Reduced density matrix 

To study the RFS, we need to know the RDM between 
two spins, and through the whole discussion we restrict to 
the case of two NN spin pairs. The Hamiltonian has the 



SU(2) symmetry, i.e. 



(7 = x,y,z) , 



which guarantees the RDM between two NN spins is of 
the form [H 



Pij = diag(gii, Q2) 



(2) 



ei = 



u^ 



92 = 



u w 



(3) 



in the basis {|00), |11), |01), |10)}, where a^\0) = 
and (Tz\l) = The matrix elements are given by [2£ 



w 



(1 ± {a^^aj^)) 



(4) 



This implies the RDM pij is only related to the spin 
correlator {(JizOjz)- It is noticed that both qi and Q2 
are Hermitian, and they can be rewritten in terms of 
Pauli operators as qi = ul, Q2 = wl + zcTx, where I 
denotes a 2 x 2 identity matrix. Therefore, it is found 
that Qi = Qi (a) (i = 1, 2) commutes with = QiiaA- 5) 
with 5 a small perturbation of the control parameter a, 
i.e., [Qi, 'qi] = 0. This commuting property will great 
facilitate our study of RFS below. 

In addition, there is an translational invariance in the 
Hamiltonian due to the periodic boundary condition, 
which leads to the fact that any two terms of the form 
(Si • Si+i) equals to each other. Applying the Feynman- 
Hellman theorem [2^, i.e., daEn — {n\daH\n) with \n) 
the non-degenerate eigenstate of Hamiltonian H and En 
the eigenenergy. To the GS, the spin correlators corre- 
sponding to two NN spin pairs are written as 



(0'lz(T2z 



:Oaeo, 



(5) 



where eg = Eq /N represents the GS energy (denoted by 
Eq) per spin. The above equation gives a direct rela- 
tion between the spin correlators and the GS energy and 
its first derivative. In other words, the elements of the 
RDMs are completely determined by cq and daeo- 



B. Reduced fidelity susceptibility 

First, we briefly review the definitions of fidelity and 
fidelity susceptibility. For two pure states |^'(q;)) and 
|4'(a -|- (5)) with d a small change of the external param- 
eter a, their overlap or fidelity is defined as [3] 



F{a) = |(*(a)|*(a + (5))|. 



(6) 



The extension to the mixed states is in general the 
Uhlmann fidelity [13, [111 

F{a) = tr^ p{ay/^p(a + 5)p{aY/\ (7) 

with p{a) and p{a -1- 5) the two density matrices. The 
fidelity susceptibility is defined as 



X = lima^o- 



-21nF 



(8) 
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Thus the fidehty susceptibihty does not depend on 6. 

Then, we will generally calculate the fidelity between 
two Hermitian and semi-positive definite density matri- 
ces g = g (a) and g = g(a + S) , which are commute with 
each other, i.e., [g, g\ = 0, so that they can be diagonal- 
ized simultaneously. With the definition of fidelity, we 
get 



F„ = trv/f?i/2^gi/2 



A,; A,; 



(9) 



where A^s and A^s are the eigenvalues of g and p, respec- 
tively. Since zero eigenvalues have no contribution to Fg , 
we only need to consider the nonzero ones. In the fol- 
lowing, the subscript i in only refers to the nonzero 
eigenvalues of g. 

For a small change S, Xi can be expanded as Xi = 
X{a + 5) ~ A, + (d^Xi) 5 + (dlX,) + O (S^^). Then 
the fidelity for matrix g becomes 



2X, 



(10) 

Here we have neglected small terms higher than second 
order. Since Y,i \ = 1, we have daXi = d^Xi = 0. 
Thus the fidelity is further reduced to 



{daX,Y 
4A,: 



(11) 



Therefore, according to the relation between fidelity and 
susceptibility, F = 1 — which is equivalent to 

Eq. ([5]), the fidelity susceptibility Xg corresponding to 
the matrix g is obtained as 



Xe 



E 



4A,; 



(12) 



This expression of fidelity susceptibility is valid for any 
commuting density matrices, and the second power on 
the right-hand side of the equation will lead to an inter- 
esting relation between the RFS and the second deriva- 
tive of GS energy shown in Eq. (fT5| . 



C. Connection to quantum phase transitions 

In the dimerized model, as the two-spin RDMs with 
different parameters commute, Eq. is applicable. By 
using the expression of the RDM (see Eqs. (O-©), after 
some calculations, the RFS for the density matrix pij is 
derived as 



(3 + 4(S,.S,))(l-4(S, -S,)) 



3 ida(cr^zcrjz)y 



4(l + (a,,a,,))(l-3(a,,a,,))' 



(13) 



which depends on both the spin correlator {aizCTjz) itself 
and its first derivative. In fact, to ensure the eigenval- 
ues of gi and g2 positive (we do not consider the zero 
eigenvalues), it is required that 



'lZ<-'jZi 



(14) 



which subsequently guarantees the susceptibility is non- 
negative. 

Now, substituting Eq. fS]) into Eq. (fT3|) . one can get 
another forms for the RFSs xi2 and X23 as follows 



X12 



X23 



16«2 (dleoY 



(3 -I- 8eo — SadaBo) (1 — 8eo -I- SadaCo) ' 



16 (dleo) 



(3 -I- 89„eo) (1 



(15) 



in terms of GS energy and its first and second derivatives. 

One key observation is that the numerators of the 
above two expressions happen to be proportional to the 
square of the second derivative of GS energy. Since the 
first derivative of energy is easily to be checked continu- 
ous (see Eq. (|21l) ') and the denominators are ensured to 
be positive and finite by Eq. (fT4|) , the singularities of the 
RFSs are determined only by the numerators. That is, if 
the second derivative of GS energy is singular at the criti- 
cal point, the RFSs is singular too. On the other hand, it 
is known that the divergence of the second derivative of 
GS energy reflects the second-order QPTs of the system, 
which is shown in Ref. [l^ explicitly as 



N 

E 



2|(^„|gi|^„) 
N{Eo - E^) 



(16) 



where Hi — dcH is the driving term of the Hamilto- 
nian iJ, and |5'„) is the eigenvector corresponding to the 
eigenvalue of H . Eq. (|16p shows that the vanish- 
ing energy gap in the thermodynamic limit can lead to 
the singularity of the the second derivative of GS energy. 
Therefore, both the two-spin RFSs can exactly reflects 
the second-order QPTs in this model. In addition, the 
second power in the numerators of the expressions, which 
origins from the relation obtained in Eq. (|12p. indicates 
that the two-spin RFSs is more effective than the second 
derivative of the GS energy in measuring QPTs. Further- 
more, by the fidelity approach, it will be shown in Sec. 
Ill that the dimerized AHC has a second-order critical 
point at a = 1. 



III. FINITE-SIZE AND CRITICAL BEHAVIOR 

In this section, we consider the critical behavior of the 
ID spin-1/2 dimerized AHC. It is known that, for < 
a <C 1, the coupling between two dimers is so weak that 
all the spins are locked into singlet states, while for a = 1, 
the system is reduced to the uniform AHC. Hence, it 
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FIG. 1: Reduced fidelity susceptibilities (a) xi2 and (b) X23 versus a for A*' = 6, 8, 10, 12 for the dimerized Heisenberg chain. 



has already been proved that the dimerized AHC has a 
critical point at a = 1 [lU, [H [H, [131, which exactly 
exists in the thermodynamic limit N oo. 



A. Finite-size behavior 

1. Analytical results for N — 4 case 

For the case that the total spins = 4, the analytical 
results can be obtained. In this case, the GS energy per 
spin of the system is [1^ 



eo 



1 / 1 + a / r\ 



(17) 



with its first and second derivatives being 



dleo = 



1 - 2a 



Vl — a + 
3 



16 (1 - a + a2) 



3/2 ■ 



(18) 



Then the susceptibilities of the RDMs pi2 and p23 can 
be derived from Eq. (fT5)) as 



(19) 



Xl2 — X23 — 



16 (1 - a + a2 



From Eq. (fT9|) we see that xi2 and X23 have the same 
expressions, and there is no singularity over parameter a. 
However, take derivation of the expression with respect 
to a, one will find that there is a maximum of X12 (or 
X23) at a = 0.5, which is also the maximum position 
of d^eo as shown in Eq. (jlSp . However, the maximum 
position a = 0.5 deviates from the real critical point 
a = 1 and can be called pseudo- critical point due to the 
finite size of the system. In addition, the different powers 
in the expressions of X12 (or X23) and O^cq over the factor 



(1 — a + a^), i.e., the former is 3/2 and the latter is 2, 
shows that the RES is more sensitive around the critical 
point. 

Besides, the exact equivalence between xi2 and X23 
is in contract with concurrences as shown in Ref. [23|. 
There, the concurrences for the reduced system, i.e., C12 
and C23 are unequal to each other and have a crossing 
point at a — 1, which leads to the mean concurrence 
takes its maximum at the critical point a = 1. This is 
because the concurrences C12 and C23 are only related to 
the GS energy and its first derivative over a, respectively. 
However, the RESs shown in Eq. ((T5)) are also determined 
by the second derivative of GS energy, which leads to the 
identical behavior between xi2 and X23- 



2. Numerical results for N — 6,S, 10, 12 



Eor the case that the total spins A^ > 4, we use exact 
diagonalization method to examine the critical behavior 
of the system in terms of the RESs. The results for A^ = 
6, 8, 10, 12 are shown in Eig. [TJ 

It is seen that both the RESs X12 and X23 can well re- 
flect the critical behavior of the system. With increasing 
system size, the pseudo-critical point exhibited by X12 
(or X23) approaches to the real critical point a = 1. Be- 
sides, the larger A^ becomes, the higher and shaper the 
peak of X12 (or X23) is- 

It should be noticed that there is a slight difference 
between xi2 and X23 for a given a and A^, which results 
from the difference between the spin correlators shown in 
Eq. In fact, the two spin correlators are equivalent, 
if we exchange the two kinds of NN couplings. Thus X12 
and X23 are also equivalent in identifying QPTs. 
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FIG. 2: Reduced fidelity susceptibifities versus a in the ther- 
modynamic limit for the dimerized Heisenberg chain, with 
c = 0.3891 and the exponent p = 1.4417 over the range of 
0.001 < 77 < 0.1 [13, i.e., 0.8182 < a < 0.9880. 



B. Infinite-size critical behavior 

Now, we consider the thermodynamic limit. To be con- 
sistent with the former works, we adopt a new parame- 
ter 77 = (1 — a)/(l-|-a). When the system approaches 
to the uniform chain limit, i.e., rj — > 0, analytical studies 
obtained by renormalization group [SO, HH had predicted 
that the GS energy per spin eg should diverge as a power 
law times a logarithmic correction, i.e., 77 "^/^/j ln?7|. How- 
ever, it is restricted to an extremely small range r] < 0.02 
[2^ . Thereafter, some numerical results pointed out that 
a pure power-law behavior is reasonably simple and ac- 
curate for larger 77 as well [H, IH, HI] . 

For generality, we assume a power law of eg as the 



form cr/P with c an overall constant. The exponent p 
are given differently over different 77 ranges. Hitherto, 
almost all the works [l^, [H, [s^] show that 1 < p < 2 
over the range < 77 < 1. For example, using the DMRG 
approach, in the exponent is fit to be p = 1.45 over 
the ran^e of 0.008 < 77 < 0.1 with c = 0.39, and in 
Ref. [33|, it is estimated in the range of 0.001 <77< 0.1 as 
that p = 1.4417 with c = 0.3891. Thus we will restrict 
1 < p < 2 in the following. The GS energy per spin in 
the thermodynamic limit can be written accordingly as 

M 

eo(^) = T^(eo(0)-c7;P), (20) 
1 + 7; 

where eo(0) = 1/4 — In 2 is the GS energy per spin for 
77 ^ 0. 

The above expression shows that the GS energy follows 
the power law behavior 77^. This gives a prediction of the 
critical point of the RFSs. From Eq. ((20|) . we can easily 
get the first and second derivatives of GS energy per spin 
in the thermodynamic limit as 

a„eo = ^{2p + a^l){l + a)-^{l~af-\ 
dleo = -2c{p-l)p{l + a)-^''+''^ {l-af-\{21) 

It is seen that, as a > and 1 < p < 2, the first derivative 
of Co does not diverge for any allowed a value, while 
the second derivative of eo has a singular point a — 1. 
According to Eq. (fT5|) . it is no doubt that the RFSs also 
diverges at a = 1. That is, the dimerized AHC has a 
second-order critical point a — 1. 

Next we discuss the critical behavior of the RFSs 
around the critical point. Insert Eq. ((2T|) into Eq. (|15p . 
we obtain the RFSs as 



Xl2 = 



X23 = 



16 



16 



[p + rj- pr]) ifP + c (2 In 2 - 1) (p + 77 - p?;) 77I+P + In 2 (In 2 - 1) 77^ 



{p-ri + prj) 772P - c (2 In 2 - 1) (p - 77 + pr]) r]^+P + In 2 (In 2 - 1) 77^ 

I 



(22) 



When a — s- 1, i.e., 77 — > 0, we only consider the leading 
terms in the expressions and get the critical behavior of 
the RFSs as 



Xl2,X23 ~ V 



2p-4 



(1 



\2p-4 



(23) 



Obviously, for 1 < p < 2, both of them diverge at 77 0, 
i.e., a = 1, as displayed in Fig. O It is shown that 
the two RFSs diverge quickly when a approaches to 1. 
For a given a, xi2 and X23 are remarkably larger than 



those in the finite-size cases. In addition, the different 
power between d^eo and Xi2(X23) over the factor (1 — a) 
indicates that these RFSs are more singular around the 
critical point. 
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IV. GENERAL CONNECTION BETWEEN 
REDUCED FIDELITY AND QUANTUM PHASE 
TRANSITIONS 

A. More exemplifications 

In the above, we have illustrated the connection be- 
tween RFS and QPTs in the spin- 1/2 dimerized AHC, 
which has SU(2) symmetry. Actually, it is straightfor- 
ward that the RFS expression (fT3|) is general for an arbi- 
trary spin-1/2 Hamiltonian with SU(2) symmetry. Com- 
bined with the translation symmetry, it is easy to obtain 
a relation between RFS and QPTs hke Eq. pS]) . In the 
following, we would like to give another two exemplifi- 
cations, which are also of SU(2) and translation symme- 
tries. One is the mixed-spin (1/2, S) dimerized Heisen- 
berg chain with S an arbitrary spin length, the other one 
is the spin-1 bilinear- biquadratic model. 

The Hamiltonian for the mixed-spin dimerized Heisen- 
berg chain with alternated spins Si and S2 is 



N/2 



Hp = ^ (Si,j-S2,j + aS2,j-Si,i+i) , 



(24) 



i=i 



here, Si and S2 denote the spin-1/2 and spin-S' operators 
respectively, and a is the ratio between the two kinds of 
NN spin couplings. The periodic boundary condition is 
assumed. As the system is of SU(2) symmetry, the RDM 
between two NN coupling spins can be expressed in the 
coupled angular momentum representation as [s^l 



F 



S-1/2 



= ^ E \s - 1/2, J.){S- 1/2, J, 

.7= = -S+l/2 

1-F 



2S + 2 



J2 \s + 1/2, j,){s + 1/2, J,\, 



J3=-S-l/2 



(25) 



with Jz the total spin along the z direction of the two 
spins and F = {S - 2(Si,j-S2j)) /(2S' + 1). Since p^j is 
diagonal, the RFS expression (|12p is available. Then 



X^j 



(a„(Si,,.S2,,))' 



(5-2(Si.,.S2j))(5 + l + 2(Si,,.S2j))' 



(26) 



Meanwhile, the system is translational invariant. Thus 
applying the Feynman-Hellman theorem to the GS of the 
system, we get the expressions for the two kinds of RFSs 
between two NN spin pairs 



X12 



X23 



(dleo) 



{S - 4eo -I- 4a9aeo) (5 + 1 + 4eo - 4a5aeo) ' 



4 (a^eo) 



{S - Ada^eo) iS + l + 49„eo) ' 



(27) 



Obviously, when S = 1/2, the above expression reduces 
to Eq. (fT5|) . The two RFSs is proportional to the second 
derivative of the GS energy per spin cq = eo(a). That 
is, the RFS also has the possibility to signal the second- 
order QPTs of a mixed-spin system. 

Furthermore, the general expression of RFS (fT3|) could 
also be extended to high-spin system, like the spin-1 
bilinear-biquadratic model, which describes the structure 
of some materials, such as LiVGe^Ofi (stI. Issjl . The Hamil- 
tonian reads 



H 



N 

BB^J2 [cos0(SrS,+i) -|-sin0(SrS,+i)2] 

i=l 



(28) 



here Si denotes the spin-1 operator, and 9 reflects the 
different coupling strengths. The periodic boundary con- 
dition is assumed as well. Obviously, this Hamiltonian is 
also of SU(2) and translation symmetries. In Eq. (24) of 
Ref. [3^ , the QPT of this model is studied by using the 
RFS between NN-coupling spins, which happens to be 
proportional to the second derivative of the GS energy 
density eo=eo(6'), i.e.. 



X12 oc (eo -I- dgeof 



(29) 



This further confirms that the two-spin RFS is an effec- 
tive tool to reveal the second-order QPTs even for high- 
spin systems. 

All the above results show that the two-spin RFS 
is tied to the corresponding spin-correlator, while the 
latter is related to the GS energy through Feynman- 
Hellman theorem. Then the RFS is connected to the 
square of the second derivative of the GS energy, which 
is used to characterize the second-order QPTs, as shown 
in Eqs. (fT5|) . ((27|l . ([29ll . Moreover, the square relation 
between the two-body RFS and the second derivative of 
GS energy Eq. (fTS]) holds for a broad class of systems 
with SU(2) and translaton symmetries, and thus the RFS 
might be more sensitive than... the second derivative of 
the GS energy in characterizing QPTs. 



B. General connection between reduced fidelity 
and quantum phase transitions 

The models considered above are all of SU(2) and 
translation symmetries. It is noticed that the definition 
of RF ([7]) depends only on the RDM, which may con- 
tains sufficient information about QPTs. This inspires 
us to infer that for a more general case, QPTs are essen- 
tially related to the RDM. In [s^, they have provided 
a powerful substantiation. They demonstrated that, un- 
der certain general conditions, the elements of two-body 
RDM are able to signal the QPTs. They consider a gen- 
eral Hamiltonian that contains two-body interaction like 



Cti) iPi 



(30) 
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where i, j enumerate N particles and is a basis for 

the Hilbert space. For the nondegenerate GS its GS 
energy is Eq — {il;\H\ip), and the element of the corre- 
sponding two-particle RDM is p!;'^^^ = {ip\aiPj){%Kj\ip). 
Thus the relation between energy and RDM is Eq ~ 

j:,MU{ij)p% where U^p.M) = e^,^^./^. + V^^j,^^ 
witn Ni the number of particles that particle i interacts 
with and the Kronecker symbols on particle j. Then 
using the Feynman-Hellman theorem, the derivatives of 
energy per particle {eo=Ea/N) are obtained as 

= j^J2^'ii^^ump^J], (31) 

ij 

+tr[{d^Umd^p,,]}, (32) 

where it follows from Eq. ([3T|) that J2ij^^i^i'''j) i'^iPij)] — 
0. As is known, according to the classical definition of 
phase transitions given in terms of the free energy ^40] , in 
the limit of T = 0, a first-order QPT (second-order QPT) 
is characterized by a discontinuity in the first (second) 
derivative of the GS energy (see also Eq. (|16p ). There- 
fore, if U{ij) is a smooth function of the Hamiltonian 
parameter ^, the origin of first-order QPTs is due to the 
discontinuity of one or more of the Pij's at the critical 
point according to Eq. ipij) . Whereas, if pij is finite at 
the critical point, the origin of second-order QPTs is the 
fact one or more of the daPij^s diverge at the critical 
point. 

Based on these facts, one find that if U{ij) is a smooth 
function and the first derivative of the elements of pij di- 
verges at the critical point, then 9|eo diverges too, which 
indicates a second-order QPT. For example, in our mod- 
els considered, the elements of pij are decided by {(TizCfjz), 
which is connected to d^e^ via Feynman-Hellman theo- 
rem. Thus the relation revealed by Eq. ([5^ may be the 
origin of the relation between RES and QPTs. This is 
not restricted to the systems with SU(2) invariance, and a 
more explicit and direct relation between RES and QPTs 
may need further deep considerations. 

In addition, the relation between the reduced fidelity 
(denoted as Fr) and its corresponding global fidelity Fq 
is given already as Fq < Fn [4l[. According to relation 
([5]), the corresponding susceptibilities satisfy yg > Xr- 

However, all the previous works [i3, [H, [il o m, 

and this work confirm that the reduced fidelity approach 
is as effective as global fidelity in characterizing QPTs, 



and in some cases, such as the models mentioned above, 
it is only necessary to know the GS energy of system in 
calculating the RES, rather than its GS for the global 
fidelity, which is generally not easy to be obtained. 



V. CONCLUSION 

In conclusion, we have studied the critical behavior 
of the ID spin-1/2 antiferromagnetic Heisenberg chain 
with dimerization in terms of RES. For the GS of the 
system, two kinds of RESs between two NN spin pairs are 
considered. It is interesting that, due to the SU(2) and 
translation symmetries, the singularities of these RESs 
are just determined by the square of the second derivative 
of the GS energy, which means the RES is more effective 
than the second derivative of the GS energy in identifying 
the second-order QPTs. Explicit calculations are also 
carried out both in finite-size and infinite-size situations. 
It is found that, as the system size increases, the pseudo- 
critical points of the RESs approach to the real critical 
point a = 1. In the thermodynamic limit, we give the 
critical exponent of the two RESs. These results further 
convince us that the critical behavior of the system can be 
reflected by the fidelity of its two-spin subsystem, which 
is of practical use in experiments. 

Furthermore, concluding the results obtained in the 
spin-1/2 dimerized model, we examine another two ex- 
amples, i.e., the mixed-spin dimerized Heisenberg chain 
and the spin-1 bilinear-biquadratic model, which have 
the SU(2) and translation symmetries as well. It is also 
found that the RESs are directly connected to the square 
of the second derivative of the GS energy, which indicates 
that for a broad class of systems with SU(2) and transla- 
tion symmetries, the RES is more effective than the sec- 
ond derivative of the GS energy in reflecting second-order 
QPTs. Moreover, the origin between RES and QPTs is 
found to be generally rooted in the relations between the 
elements of the RDM and derivatives of GS energy, which 
provides a general proof for the fact that RES can be used 
to detect the second-order QPTs of the system. 
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